Msat = 900000
arm_Separation := 76e-9
nm_per_turn_d := 350
Aex = 1e-11
alpha = 0.5
MinimizerStop = 1e-1
MinimizerSamples = 20
FixDt = 1e-13
Mesh_size := 0.0
tableAdd(E_exch)
tableAdd(E_demag)
tableAdd(E_total)
tableAdd(LLtorque)
tableAdd(MaxAngle)
tableAdd(Aex)
TableAddVar(Mesh_size, "Mesh_size", "m")
tableAdd(Edens_exch)
tableAdd(Edens_demag)
tableAdd(Edens_total)
tableadd(B_ext)
OutputFormat = OVF2_Binary
edgesmooth = 8
SetMesh(216/2,216/2,1250/2, 2e-9,2e-9,2e-9, 0, 0, 0)
m.loadfile("~/Double_Helix_Hybrid_Vortex-AP_State_MinorLoop.ovf")
save(m)
tablesave()
Bmax := -70.0e-3
Bstep := 0.5e-3
MinimizerStop = 1e-6
ang := 15 * pi / 180
for B := 27.0e-3; B >= Bmax; B -= Bstep {
	B_ext = vector(B*sin(ang), 0, B*cos(ang))
	relax()
	tablesave()
	save(m)
}
